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Several  control  function  Int e.-po  1  at  1  or. 
cr.r.iqjes  in  a  general  three-dimensional  e..iptlc 
d  generation  code  and  tnelr  effects  on  f.o. 
utidr.s  using  an  implicit  Euler  algorithm  are 
amlnet.  Tr.ese  results  will  serve  to  guite  tr.e 
sign  cf  control  function  procedures  a no 
terpoiation  tec.nr.iques  in  general  grid 
neratlon  codes.  Tnree  configurations  and  t.oree 
id  types  !"C",  “C",  and  "H"  grids  ,  are  examined. 
e  results  Indicate  tnat  tne  selection  of  ire 
r.trcl  function  Interpolation  techniques,  wr.lcr. 
fects  grid  spacing,  should  be  based  on  boundary 
■vature  and  spacing.  Tne  selection  of  tne 
terpolation  tecnnlque  can  then  be  made 
ar.sparent  to  the  user  of  general  grid  generation 


Nomenclature 

elements  of  covariant  metric  tensor 
tm  =  1,2.3;  n  =  1,2,3) 
elements  of  contravarlant  metric  tensor 
U  =  1,2,3!  n  =  1,2,3) 

Cartesian  position  vector 

(r  r  xl  *  yj  .  zk) 

three  Cartesian  coordinates 

three  Curvilinear  coordinates  (1  = 

1,2,3)  of  transformed  region 

"Control  functions"  which  serve  to 

control  spacing  and  orientation  of  grid 

lines  (n  =  1,2,3) 

local  radius  of  curvature 

diameter 

freestream  Mach  number 
angle  cf  attack  (degrees) 
pressure  coefficient 
sonic  pressure  coefficient 

Introduction 

e  field  cf  computational  fluid  dynamics 
•rent  research  Is  aimed  primarily  at 
both  the  accuracy  and  efficiency  of  the 
techniques  employed  in  solving  fluid 
i  transfer,  and  combustion  problems.1 
r.ts  m.st  be  made  in  both  of  these  areas 
tc  become  a  viable  partner  to  tne  wind 
aerody  >  u  !  c  design.  Numerical  grid 
r.  requires  tecnnoiogy  development  and  has 
a  as  a  major  pacing  Item  for  realistic 
missile  applications.  Tne  techniques  of 
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numerical  grid  generation,  and  Its  application  In 
tne  numerical  solution  of  pa-tiai  differential 
equations,  are  covered  Ir.  detail  in  a  recent  text 
or.  tne  subject.-^  Several  surveys  of  the  field 
nave  also  been  given, 1  and  three  conference 
proceedings  dedicated  to  this  area  have  been 
putlisned.0-'1' 

Tne  computational  grid,  or  mesh,  working  In 
conjunction  with  the  flow  solver  can  have  a 
substantial  Impact  on  the  aerodynamic  solution;1 
therefore,  It  Is  the  objective  of  this  analysis  to 
examine  several  grlddlng  procedures  and  tnelr 
effects  on  Euler  approximations.  This  study  was 
conducted  as  part  of  the  development  effort  of  a 
general  three-dimensional  grid  code  for  realistic 
alrcraf t/mlssl le  configurations.  The  purpose  of 
tnls  analysis  Is  to  demonstrate  the  viability  of 
automatic  determination  of  the  control  function, 
in  the  elliptic  grid  generation  system,  from  the 
boundary  point  distributions  using  separate 
interpolation  techniques  for  the  terms  arising 
from  spaclhg  and  local  curvature.  These  results 
will  serve  to  guide  the  design  of  control  function 
procedures  and  interpolation  techniques  to  be  made 
automatic  and  transparent  to  the  user  in  this  code 
and  are  similarly  of  relevance  to  the  design  of 
general  grid  generation  codes  by  others. 

Grid-Generation 

A  general  three-dimensional  grid  generation 
code  has  been  written  which  allows  any  number  of 
blocks  to  be  used  to  cover  an  arbitrary  three- 
dimensional  region.  Any  block  can  be  linked  to 
any  other  block,  with  complete  continuity  across 
the  block  interfaces.  The  composite  structure  is 
such  that  completely  general  configurations  may  be 
treated,  the  arrangement  of  the  subregions  being 
specified  by  Input,  without  modification  of  the 
code.  Tne  code  Includes  a  three-dimensional 
algebraic  generation  system  based  on  transflnite 
Interpolation  for  the  generation  of  an  Initial 
solution  to  start  the  Iterative  solution  of  the 
elliptic  generation  system.  This  feature  also 
allows  the  code  to  be  run  as  an  algebraic 
generation  system  If  desired. 

Tnls  code  uses  an  elliptic  generation  system 
with  automatic  evaluation  of  control  functions 
from  boundary  point  distributions.  Several 
procedures  for  determining  the  control  T—nct  lens 
from  the  boundary  point  distributions  are 
Incorporated  in  the  oode,  and  it  was  the  purpose 
of  the  present  study  to  evaluate  the  effectiveness 
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of  tneae  procedures  with  regard  to  Euler 
solutions . 
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Trie  elliptic  grid  generation  system  is  defined 

by: 
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wn ere  the  g31”  are  the  elements  of  the 
contra  variant  metric  tensor: 
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so  that 

p  .  V  (*t* 

i  ~~v~  (6) 
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With  the  J  points  on  the  vertical  sides,  and  I 
pcints  on  the  horizontal  sides,  the  control 
function  on  the  vertical  sides  then  can  be 
evaluated  musing  a  central  difference  technique 
for  Eq.  (6)  )  from  the  point  distribution  thereon 
as 

P  ( 1  ] )  -  -  I  r  ( 1 .  J  *  1 )  -  z  T  ( 1 .  i )  *  y  ( 1 .  J  -  i ) : 


These  element:  are  more  conveniently  expressed  In 
terms  of  the  elements  of  the  covariant  metric 
tensor,  s_„: 


S--  =  r  (  i  •  _r  f  n  *  -  (3) 

wrier.  can  be  calculated  directly.  Thus 

l**  *  K*,k  -«,!«;«>  '«) 

i 

(m,i,J)  cyclic,  (n,k,l)  cyctlo, 
where  g,  the  square  of  the  Jacobian,  Is  giver,  by 


8  =  del  18,,  !  =  -V  '  (t(l  x  -r(3  1  (5) 

Ir.  t.nese  relations,  _r  is  the  Cartesian 
position  vector  of  a  grid  point  (r  s  lx  ♦  Jy  . 
kz),  and  t he  «*  (1=1, 2,3)  are  the  three 
transformed,  curvilinear  coordinates.  The  Pn  are 
the  "control  functions"  which  serve  to  control  the 
spacing  and  orientation  of  the  grid  lines  In  the 
field. 

Negative  values  of  the  control  function  Pn 
cause  grid  lines  on  which  (n  is  constant  to  move 
In  the  direction  of  decreasing  (n.  This  feature 
can  be  used  to  concentrate  grid  lines  near  other 
grid  lines  and/or  points  or  In  certain  regions  of 
physical  space.  However,  a  more  automatic 
procedure  Is  to  dete:mlne  the  oontrcl  functions  so 
as  to  project  the  boundary  point  spacing  Into  the 
field.  The  details  are  discussed  In  reference  3; 
however,  a  brief  explanation  Is  presented  here. 

Consider  first  a  rectangle  with  equally  spaced 
points  on  the  horizontal  sides  but  unequal  spacing 
on  the  two  vertical  sides  (Fig.  1).  With  no 
control  functions,  i.e.,  Pn  =  0,  Eq.  (1)  will 
produce  a  grid  that  attempts  to  be  equally  spaced 
In  the  interior  of  the  region,  (Fig.  2).  A  grid 
of  parallel  lines  for  this  configuration, 
reflecting  the  unequal  spacing  on  the  boundaries 
(Fig.  3),  can  qnly  be  produced  from  Eq.  (1)  witn 
(taxing  i  '  to  vary  on  tne  horizontal  sides 
and  to  vary  on  the  vertical  sides). 

Tne  proper  values  of  P2  needed  to  accomplish 
this  are  determined  by  evaluating  Eq.  (1>  or.e- 
dimenslonal  ly  cn  the  vertical  sides,  with  tr.e 
result 
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ylrU-i  +  i;  -  y(i.j-t)] 

for  J=2,3,--,J-1,  with  an  analogous  equation  with 
I  for  the  first  argument.  The  values  of  P2  In  tne 
Interior  of  the  region  then  can  be  interpolated 
from  values  of  P^  on  the  two  vertical  sides.  A 
similar  evaluation  of  the  other  control  function, 
Pj,  on  the  two  horizontal  sides  from 


produces  a  zero  In  the  present  case  due  to  equal 
spacing  on  these  sides.  In  the  case  of  unequal 
spacing,  the  values  of  P,  In  the  Interior  of  the 
region  would  be  evaluated  by  Interpolation  between 
the  values  on  the  two  horizontal  sides.  With  the 
control  functions  evaluated  In  this  manner,  Eq. 

(1)  will  produce  a  grid  composed  of  parallel 
straight  lines  for  this  boundary  configuration, 
thus  reflecting  the  boundary  point  spacing  into 
the  field  (Fig.  3). 

Now  consider  an  O-type  grid  with  two 
concentric  circular  boundaries  and  equally  spaced 
points  around  the  circles  (Fig.  4).  Because  of 
Its  inherent  tendency  to  cause  the  grid  lines  to 
move  closer  to  convex  boundaries,  Eq.  (1)  with  no 
control  functions  will  produce  a  grid  with  unequal 
radial  spacing  of  the  circumferential  lines  (Fig. 
5).  In  this  case,  evaluation  of  Eq.  (i)  with  the 
polar  coordinate  transformation; 

x ( r . «)  =  r((Z (cost 

J  !  t  ■  8)  *  rtf*)  sinP 

yields  the  equation 

Fj  =  -  ♦  r,z 


Tnus,  In  order  to  produce  a  specified  radial 
distribution  of  lines,  the  control  function  P2 
must  be  evaluate^  from  Eq.  (8)  using  the  given 
distribution  r(i‘).  First,  a  radial  out  is  made 
in  the  physical  plane  from  the  inner  to  the  outer 
circle.  The  circular  region  can  now  be  unwrapped 
to  form  a  computational  field  that  has  Its  two 
vertical  sides  corresponding  to  the  cut  between 
tr.e  two  circular  boundaries  In  the  physical  field 
\Flg.  6).  Therefore,  with  the  speolfled  radial 
distribution  placed  on  the  two  vertloal  sides  of 
tne  computational  field,  the  oontrol  function  P2 
can  be  evaluated  from  Eq.  (6)  on  theae  sloes  as 
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ana  a  similar  equation  with  I  as  the  first 
argument.  Tne  interior  values  then  can  be 
determined  by  interpolation  between  the  two 
vertical  sides  as  before.  Note  that  this  amounts 
to  interpolation  in  the  circumferential  direction 
In  tne  physical  plane.  Again,  equal  spacing 
around  the  circles  produces  a  zero  value  of  the 
otner  control  function,  ?..  In  general,  P,  would 
be  evaluated  on  the  two  circles  and  Interpolated 
between  the  two  horizontal  sides  In  the 
computational  region,  i.e.,  between  the  two 
circles  In  tne  physical  field.  Both  tne  radial 
and  ci rcumferent lal  Interpolation  schemes  are 
represented  In  Figure  7. 

Tne  second  term  In  Eq.  (8)  arises  from  the 
curvature  of  the  boundary,  and  the  denominator  Is 
tne  local  radius  of  curvature  of  the  grid  line 
that  la  to  pass  through  the  point  where  the 
control  function  Is  being  evaluated.  Tnls  term 
acts  to  reduce  the  magnitude  of  the  control 
function  in  order  to  allow  for  the  natural 
tendency  of  the  grid  lines  to  move  toward  convex 
boundaries.  Since  the  lines  tend  to  concentrate 
near  the  Inner  circle  even  with  zero  control 
functions,  the  use  of  the  first  term  alone  In  Eq. 
(8),  (In  analogy  with  the  flat  boundary  case, 
i.e.,  Eq.  (6)),  would  produce  a  stronger 
concentration  of  lines  near  the  Inner  circle  than 
was  intended. 

Finally,  consider  a  C-type  grid  (Fig.  8),  with 
the  resulting  computational  region  (Fig.  9).  If 
now  the  control  function  P,  Is  evaluated  on  lines 
1-2  and  **-3  In  the  physical  field  from  Eq.  (6), 
and  the  Interior  values  are  Interpolated  between 
the  two  vertical  sides  In  the  computational  field, 
the  resulting  control  function,  while  serving  well 
over  the  right  portion  of  the  physical  field,  will 
be  too  strong  over  the  Inner  body  where  the  line 
curvature  Is  not  zero.  The  use  of  Eq.  (8)  on  the 
lines  1-2  and  ^-3  would  be  no  better  since  the  r 
Ir.  the  denominator  Is  to  be  Interpreted  as  the 
local  radius  of  curvature  of  the  crossing  line  and 
hence  is  Infinite  on  these  lines  so  that  the 
second  term  in  Eq.  (81  vanishes. 

This  situation  can  be  remedied  by 
Interpolating  for  the  local  radius  of  curvature  in 
Eq.  (8)  between  the  Inner  and  outer  boundaries  In 
toe  physical  field,  I.e.,  between  the  horizontal 
sides  In  the  computational  region.  However,  since 
tne  .  derivatives  In  Eq.  iS)  oust  still  be 
evaluated  on  tne  vertical  sides  it  is  necessary  to 
separate  Eq.  ,9  into  three  pieces: 
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(10) 


Now  the  two  quantities  in  parentheses  are 

evaluated  on  t  r,e  vertical  sides  of  the 
:-rz;  at  at  ior.a  1  region,  while  the  radius  cf 
curvature  is  evaluated  on  the  two  horizontal 
sides,  "he  control  function  In  tne  Interior  then 


-v.  -  i  not 
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Is  evaluated  by  Interpolating  the  quantities  in 
parentheses  between  the  vertical  sides, 
Interpolating  tne  radius  of  curvature  between  tr.e 
horizontal  sides,  and  then  evaluating  P^  Pr°-  8q- 
(10)  using  these  Interpolated  values.  Note  that 
this  procedure  supplies  a  finite  radius  of 
curvature  over  the  Inner  body,  thus  reducing  the 
control  function  appropriately  In  this  region. 

A  problem  arises,  however,  when  the  radius  of 
curvature  (r)  Is  of  opposite  sign  on  the  two 
boundaries  between  which  It  Is  Interpolated  (Fig. 
10),  since  then  the  Interpolation  -ill  produce  a 
zero  value  at  some  point  In  between,  and  at  such  a 
point  the  second  term  of  Eq.  (10)  is  Infinite. 

This  problem  can  be  corrected  by  interpolating 
1/r,  Instead  of  r.  This  amounts  to  writing  Eq. 

(  10)  as 


/-r  <!<*  )  ♦  (  r ;  ± 
'  f,  2  ’  r 


(11) 


and  Interpolating  the  curvature  (1/r),  rather  than 
the  radius  of  curvature  (r). 

Although  the  exact  equations  for  the  general 
ease  are  more  complicated,  the  control  function, 

Pn,  may  still  consist  of  spacing  terms  along  a 
boundary  and  terms  arising  from  the  local 
curvature  of  crossing  lines.  This  amounts  to 
Interpolating  the  spacing  terms  between  the  f our 
sldes  on  which  fn  varies  (Fig.  11)  and 
interpolating  the  local  curvature  terns  between 
the  two  sides  on  which  Pn  Is  constant  (Fig.  12). 
The  question  that  then  arises  Is  whether  the 
transflnlte  Interpolation  for  the  spacing  term  and 
the  local  curvature  term  should  use  linear 
blending  functions  or  blending  functions  based  on 
physical  arc  length.  The  former  (linear)  amounts 
to  Interpolating  In  terms  of  the  transformed 
curvilinear  coordinate,  while  the  latter  (arc) 
amounts  to  interpolation  with  respect  to  the 
physical  distance.  For  example,  on  the  grid 
Illustrated  In  Fig.  13,  Interpolation  with 
"linear"  blending  functions  would  produce  a  value 
on  line  3  that  Is  the  average  of  that  on  lines  1 
and  5,  while  with  "arc"  blending  functions  the 
value  produced  on  line  3  would  be  closer  to  this 
on  line  1. 

Flow  Solver 

For  this  Investigation  a  three-dimensional 
time-dependent  Euler  approximation,  developed  by 
Whitfield,  ‘was  employed  using  finite  volume 
discretization  and  a  second-order  Implicit  scheme 
to  solve  the  f lux-vector-spllt  form  of  the 
equations  with  local  time  stepping.  This  scheme 
solves  5x5  block  bl-dlagonal  systems  of  equations 
using  Doolittle's  method.1^  Characteristic 
■arlable  boundary  conditions  are  used  in  the 
farfleld  and  at  impermeable  surfaces.  No 
additional  artificial  dissipation  is  added,  and 
tne  scheme  Is  conditional  ly  stable  in  three 
dimensions.  Optimum  CFL  numbers  appear  to  be 
between  12.5  to  15. 

Configuration  and  Flow  Conditions 

Three  generic  weapcn/store  configure! ions  are 
examined  using  tne  various  gridiir.g  prcoec-res  ar.d 
grid  types.  Tr.e  first  geometry  ^Configuration  ') 


3 


TT 


consists  of  a  10/3  D  cylindrical  centerbody  and 
5/3  0  tangent  ogive  forebody  and  afterbody1""  (Fig. 
1  *» )  and  is  modeled  using  an  81x15x10,  0-type  grid. 
This  configuration  Is  used  to  represent  a  typical 
fuel  tank  carried  on  current  aircraft  and  Is 
examined  at  two  transonic  Mach  numbers,  0.85  and 
1.80,  both  at  zero  degrees  angie-of-attack.  Tr.e 
second  geometry  (Configuration  2)  consists  of  a 
10°  cone-cylinder1^  (Fig.  15)  and  Is  modeled  using 
an  81x15x10,  H-type  grid.  This  geometry  Is  used 
to  represent  a  generic  penetrator  and  Is  examined 
at  ^  =  0.999  and  <r  =0°.  Th"  third  geometry 
(Configuration  3)  consists  of  a  1.5  D  tangent 
ogive  forebody  and  a  20  D  cylindrical  aft 
section  ^  (Fig.  16)  and  Is  modeled  by  a  65x20x10, 
C-type  grid.  This  conf iguratlon  Is  used  to  model 
a  hlgh-flneness  ratio  (21.5:1)  body-alone  missile 
airframe,  and  Is  analyzed  In  the  subsonic  range  at 
M^.sO.70  at  n  =5.07°  and  In  the  transonic  range  at  M 
=  0.90  at  o:0°. 

Analysis 

Two  Investigations  examined  control  function 
Interpolation  techniques.  Tne  first  analysis 
Involved  employing  both  the  "linear”  and  "arc" 
tecnnlques  for  Interpolating  the  spacing  term  In 
the  control  function.  The  second  study 
Investigates  the  use  of  the  "linear"  and  "arc" 
techniques  with  respect  to  the  radius  of  curvature 
(r)  and  curvature  (1/r)  Interpolations  for  the 
local  curvature  term  In  the  control  function. 

Spacing  Term 

Figure  17  shows  the  effects  of  the 
Interpolation  for  the  spacing  term  using  the 
"linear"  and  “arc"  techniques  for  the  C-type  grid 
on  Configuration  3.  Using  the  "linear" 
Interpolation  (Fig.  17a),  the  grid  lines  better 
reflect  the  concentrated  point  distribution  in  the 
shoulder  region  of  the  body,  thereby  aligning  the 
grid  with  the  shock  pattern  that  develops  near 
that  location.  Figure  17b  shows  a  close-up  view 
of  the  grid  near  the  shoulder  region  of  the  body. 
Using  the  "arc"  interpolation,  the  grid  lines  are 
slanted  far  forwnrd  and  hence,  do  not  align 
themselves  with  the  shock  pattern  (Figs.  17c ,d). 
Figure  18a  shows  the  effects  these  two  grids  have 
on  the  Euler  solution.  Figure  18b  gives  an 
expanded  view  of  the  shock  region.  For  the 
"linear"  grid,  a  typical  second-order  solution  to 
the  Euler  equations  at  M^.  =0.90  and  o  =0°  Is 
obtained.  One  expects  the  computed  shock  strength 
to  be  slightly  higher  than  the  experimental  data 
at  the  peak  of  the  expansion  and,  although  not 
desirable,  second-order  "ringing"  (dispersive 
effects)  Is  expected  downwind  of  the  shook  for 
this  upwind  scheme.  The  "linear"  grid  solution 
shown  has  expanded  at  the  shock  more  than  one 
expects,  but  does  reflect  the  characteristics  of 
the  invlscid  solver.  The  "arc"  grid  solution, 
although  core  closely  matching  the  experimental 
data,  does  not  reflect  the  qualities  of  an 
accurate,  second-order,  invlscid  solution.  Figure 
19  snows  the  effects  these  two  grids  have  on  the 
Euler  solution  at  M*.  =0.70  and  „  =5.07°.  The 
"linear"  grid  solution,  once  again,  reflects  a 
typical  second-order,  Invlscid  solution  by  over 
predicting  the  expansion  and,  slnoe  no  shock  was 
formed,  exhibits  no  dispersive  effeots  (ringing), 
Tne  "arc"  grid  solution  shows  similar  results  In 


which  the  expansion  Is  overpredicted  and  no 
dispersive  effects  are  witnessed  downwind  of  the 
expansion. 

Although  not  shown,  similar  results  are 
obtained  for  Configuration  1  (0-type  grid)  at  both 
flow  conditions  and  Configuration  2  (H-type  grid) 
at  =0.999  and  n  =0°.  These  results  again  show 
that  the  "linear"  Interpolation  tecnnlque  for  the 
spacing  term  reflects  a  more  accurate  second- 
order,  Invlscid  solution  for  all  cases. 

Local  Curvature  Term 

The  first  part  of  this  analysis  examines  the 
effects  of  Interpolating  for  the  radius  of 
curvature  (r)  using  both  the  "linear"  and  "arc" 
techniques  to  obtain  the  local  curvature  term  In 
the  control  function.  Figure  20  shows  the  effects 
of  these  two  techniques  for  the  0-type  grid  on 
Conf  Iguratlon  1.  Using  the  "arc"  technique,  a 
better  concentration  of  grid  lines  Is  obtained 
close  to  the  body  (Fig.  20a, b).  Using  the 
"linear"  technique,  a  comparable  grid  Is  produced; 
however,  the  grid  is  not  as  concentrated  near  the 
body  (Fig,  20c, d).  Figure  21  shows  the  effects 
these  two  grids  have  on  the  Euler  solution  at 
=  1.20  and  a  =0°.  The  "arc"  grid,  with  Its  greater 
concentration  of  lines  near  the  body,  resolves  the 
flowfield  more  accurately  In  the  stagnation 
region.  Physically  reasonable  results  are 
obtained  using  the  "arc"  grid  In  which  the  first 
cell  on  the  nose  yields  a  higher  pressure 
coefficient  (C„)  than  does  the  second  or  third 
cell  (Fig.  21b;.  The  "linear"  grid  yields 
unreasonable  results  In  which  the  third  cell  back 
on  the  nose  yields  a  higher  C_  than  does  the  first 
or  second  cell  (Fig.  21b).  Figure  22  shows  the 
effects  these  two  grids  have  on  the  Euler  solution 
at  Kao =0.85  and  «  =0°.  Both  grids  yield  physically 
reasonable,  although  markedly  different,  results 
at  these  flow  conditions. 

Similar  results  are  obtained  for  Configuration 
2  at  Me  =0.999  and  a  =0°  and  for  Configuration  3  at 
both  flow  conditions.  These  results  show  that 
interpolating  for  the  radius  of  curvature  (r) 
using  the  "arc"  technique  to  obtain  the  curvature 
term  in  the  control  function  yields  a  more 
physically  reasonable  solution  for  all  cases. 

Due  to  the  fact  that  complications  can  arise 
when  Interpolating  for  the  radius  of  curvature 
(r),  a  seoond  examination  Is  performed  to  study 
the  effects  of  interpolating  for  the  curvature 
(1/r)  using  both  the  "linear"  and  "arc" 
techniques.  Figure  23  shows  the  effects  of  using 
the  "arc"  technique  on  the  0-type  and  C-type 
grids.  This  technique  yields  an  undesirable  grid 
in  which  the  lines  are  highly  skewed  and,  hence, 

Is  not  useful  to  the  Euler  solver.  Therefore, 
when  Interpolating  for  the  curvature  (1/r)  only 
the  "linear"  technique  should  be  used. 

Control  Function 

Depending  on  the  geometry,  there  are  two 
optimum  approaches  for  obtaining  the  control 
function,  Pn.  The  first  method  1»  to  employ  the 
"linear"  teonnlque  to  obtain  the  spacing  term  and, 
when  the  two  opposing  boundaries  have  ourvsture  of 
the  same  sign,  to  Interpolate  for  the  radius  of 


curvature  (r)  using  the  "arc"  technique  to  obtain 
the  local  curvature  term  (Method  1).  The  second 
method  Is  to,  again,  employ  the  "linear"  technique 
to  obtain  the  spacing  term  and,  when  the  two 
opposing  boundaries  have  curvature  of  the  opposite 
sign,  to  Interpolate  for  the  curvature  using  the 
"linear"  technique  to  obtain  the  local  curvature 
term  (Method  2). 

These  two  methods  yield  similar  grids  and  a 
thorough  comparison  shows  that  both  methods  give 
approximately  the  same  Euler  solutions.  Figure  2M 
shows  the  effects  of  these  two  methods  on 
Configuration  1  (0-type  grid)  at  M^l.20  and  „  =  0°. 
Tne  results  are  essentially  the  same  except  in  the 
stagnation  region  where  Method  2  (spacing: 
"linear”,  local  curvature:  curvature  -  "linear") 
does  not  concentrate  the  grid  lines  as  close  to 
the  body  and;  nence,  does  not  appropriately 
resolve  the  flowfleld  (Fig.  2Mb).  Figure  25  snows 
the  effects  of  these  two  methods  on  Configuration 
2  (K-type  grid)  at  M^.  =0.999  at  a  =0°.  These 
results  show  that  both  methods  yield  essentially 
Identical  solutions.  Figure  26  shows  the  effects 
of  these  two  methods  on  Configuration  3  (C-type 
grid)  at  M^  =0.90  and  „  =0°.  Again,  these  results 
show  that  both  methods  yield  the  same  second- 
order,  lnvlscld  solution. 

As  can  be  seen  by  a  comparison  of  Figures  20 
and  23  and  the  results  In  Figure  2M,  it  is 
advisable  to  Interpolate  for  the  radius  of 
curvature  (r),  rather  than  the  curvature  (1/r),  to 
obtain  the  local  curvature  term  In  the  control 
function  when  the  curvature  has  the  same  sign  on 
the  two  opposing  boundaries.  In  fact, 
Interpolating  for  the  curvature  (1/r)  does  not 
produce  the  exact  for-m  given  by  Eq.  (8)  for  the 
simple  case  of  two  concentric  circles. 

Conclusions  and  Recommendations 

This  study  has  demonstrated  the  viability  of 
automatic  determination  of  the  control  function 
from  the  boundary  point  distributions  using 
separate  Interpolation  of  the  terms  arising  from 
spacing  and  those  arising  from  local  curvature. 

Tne  following  Is  recommended  as  the  optimum 
approach  for  obtaining  the  control  function  In  the 
elliptic  grid  generation  system  (Eq.  1).  For  the 
spacing  term,  the  "linear"  blending  approach 
should  be  employed  globally.  For  the  local 
curvature  term,  global  Interpolation  for  the 
radius  of  curvature  (r)  using  the  "arc"  technique 
should  be  employed.  However,  when  the  two 
opposing  boundaries  have  curvature  of  the  opposite 
sign,  local  interpolation  for  the  curvature  (1/r) 
using  tne  "linear”  tec.nnlque  should  be  used. 

It  has  been  shown  throughout  this  analysis 
tr.at  various  grid  generation  techniques  employed 
within  the  code  do  have  a  substantial  Impact  on 
tne  solution  and,  when  care  is  taken  to  generate 
tne  grid  In  the  proper  manner,  these  grids  help 
tne  flow  solver  to  yield  accurate  second-order, 
Ir.vlsoid  solutions. 
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Figure  1.  Unequal  Spacing  on  Vertical  Side* 
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Figure  3.  Grid  of  Parallel  Line* 


Figure  4.  F.qually  Spaced  Point!  Around  Circle* 


Figure  S.  Unequal  Radial  Spacing  of 
Circumferential  Line* 
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Figure  7.  Radial  (P^  )  and  Circumferential  (!’ 
Interpolation  Scheme* 
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Figure  8.  C-Type  Grid.  Physical  Region 


Outer  Itoundarj 


Figure  12.  Two  Sides  on  Which  P  is  Constant 
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Figure  9.  Computational  Region 


Figure  10.  Radius  of  Curvature  (r)  With 
Opposite  Signs 


Figure  13.  Interpolation  Techniques 
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Figure  14.  Configuration  1 


Figure  11.  Four  Sides  on  Which  t  Varies 


Figure  15.  Configuration  2 


Figure  16.  Configuration  3 
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Figure  17*.  'Linear'  Interpolation  for 
Spacing  Term  (C-Type) 
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Figure  17b.  Clo»e  Up  View  of  Shoulder 
Regioa  ('Linear*) 
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Figure  17c.  'Arc*  Interpolation  of 
Spacing  Term  (C-Type) 
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f  igure  17d.  Clo»e  Up  View  of  Shoulder 
Region  ('Arc') 


<  i  m  in  '  n 


Figure  18*.  Surface  Preuure  Diitribution 
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Figure  18b.  Expanded  View  of  Shock  Region 
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Figure  19.  Surface  Pressure  Distribution 


Figure  20*.  'Arc'  InterpoUtion  for  Local 
Curvature  Term  (O-Type) 


f  igure  20b.  Close  Up  View  of  Stagnation 
Region  ('Arc') 


Figure  20c.  "Linear*  Interpolation  of  Local 
Curvature  Term  (O-Type) 


Figure  20d.  Close  — Up  View  of  Stagnation 
Region  ('Linear') 
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I  igurc  21a.  Surface  Pressure  Distribution 
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Figure  24b.  Expanded  View  of  Stagnation  Region 
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Figure  26.  Surface  Pressure  Distribution 


PPCSiUrl  [ISTFiBUTlO't 
‘C,,:  IGuPPT !QN  t  ulo^e -Cy  I  mder  i 


- j - 1 

■ 

1 

■  /  \ 

«* 

X  1 

« 

• 

v. 

•  - . - 

4 

»  1 

,  £  *  •  *  ^  j~x  ^  x 

a1  ic 

!  *  ! 

1* 

iTt-’.  [ 

~b  ■'rt.f' ,  •  t?  iri* 

-o.*  2  -m*  ■  e  *1 

j 

><  d  1 

m  s  s  '■+■*  "nr.'oc 

*  :  .i  "  nr-  •  o|  CbiJ  1  ■+  '•  J: 

•  *v 

f  **v  i  f  ** 

.*  L1i» 

Figure  25.  Surface  Pressure  Distribution 


